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Abstract 

Consider a string of N + 1 damped oscillators moving in R of which the 
motion of the first (called the "leader" ) is independent of the others. Each of 
the followers 'observes' the relative velocity and position of only its nearest 
neighbors. Inasmuch as these are different from 0, this information is then 
used to determine its own acceleration. Fix all parameters except the number 
N in such a way that the system is asymptotically stable. Now as N tends 
tends we consider the following problem. At t = the leader gets kicked 
and starts moving with unit velocity away from the flock. Due to asymptotic 
stability the followers will eventually fall in behind the leader and travel 
each at its own predetermined distance from the leader. In this note we 
conjecture that before equilibrium ensues, the perturbations to the orbit of 
the last oscillator grow exponentially in N except when there is a symmetry 
in the interactions and the growth is then linear in N. There are two cases. 
We prove the conjecture in one case, and give a strong heuristic argument in 
the other. 
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1. Introduction 

In this note we study how a long, but finite, string of asymmetrically 
coupled damped oscillators reacts as one of its members (the leader) changes 
its velocity suddenly. The aim is to study how the stability of large flocks 
depends on the kind of interaction between the individual agents. 
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Figure 1: The communication graph of the system. Each agent is linearly coupled to its 
nearest neighbor. The arrows give the direction of the information flow. The interactions 
for the boundary agents are different from those in the interior. At t = the agent labeled 
0, the leader, undergoes a forced motion: a kick in the direction of the arrow above it. It 
receives no feedback from the flock. 

Assume then that we have a string of oscillators as depicted in Figure 
[U where in principle each agent observes the relative distance and velocity 
of its immediate neighbors and uses those to compute its own acceleration. 
The types of interaction are indexed by a parameter p e [0,1]: each agent 
multiplies the information coming from the neighbor 'in front' by weight 1 — p 
and from the neighbor 'behind' by the weight p. Naturally the interactions of 
the first and last agents are a little different from those with two neighbors. 
The leader itself is assumed not at all influenced by the rest of of the flock 
(an 'independent leader'). At t = the agent labeled 0, the leader, under- 
goes a forced motion: namely its velocity goes instantaneously from to 1 
(it receives a 'kick'). The response of the last agent is called the 'impulse 
response' function. We study this response function holding all parameters 
fixed, except the number N of agents. 

The interplay between graph theory and coupled linear ODE's is described 
in detail in [l|. In particular, the collection of agents connected by the 
arrows that give the information flow gives a directed graph, known as the 
communication graph, see Figured] 

In @,i we studied what happens if the weights are equal, that is: p = n- 
In that case the perturbation propagates from the leader throughout the 
flock and is roughly multiplied by N, the size of the flock, when it arrives at 
the trailing agent (labeled N). The perturbation then decays as the system 
is assumed to be asymptotically stable. In 0] and this paper we take up 
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the study of this system when p ^ |. In the former we concluded that the 
frequency response function grows exponentially for all p € [0,1). Here we 
we show that if p G (§,1), then the impulse response function also shows 
exponential growth in N. When p e (0, \) we conjecture that this is also 
true. (When p = |, both grow only linearly with N). 
We study systems of the form (for details see [3]): 



Mz + r (t) 



(i) 



where z = (z±, zi,...,Zn, ijy) is the vector of position and velocity of the 
agents. The leading car is not encoded since its orbit Zo(t) is a priori given. 
The matrix M is defined in terms of the Kronecker product (<E>) 



M = I <g> A+ P <g> K 



(2) 



Here I and P are TV-dimensional square matrices, where / is the identity and 
P is given by 



( 

1-p p 



P = I — Qp where Q, 



P 



\ 



1-p p 
1 



(3) 



The 2x2 matrices A and K are given by: 
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Finally: 
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/ 



and K 





/ g 



(4) 



o \ 

P) (fz (t)+gzo(t)) 




(5) 



V / 

We consider the problem where the flock is at equilibrium for time t < 0, 
that is: the agents are at rest and properly spaced. At time we provide 
an impulse to the leader: i' (t) = So(t), the leader then advances at constant 
velocity equal to 1 for time t > 0. This problem is ultimately motivated by 
what is called the canonical traffic problem in 0,0], in which one imagines 
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a long row of cars waiting for the traffic light to turn green. When that 
happens, the first car, the leader, quickly accelerates to the desired speed 
and the others aim to follow it. The problem now is to deduce the motion of 
the other agents Zk(t), k = 1, N and t > 0. 

The control parameters f,g<0, are assumed negative, implying asymp- 
totic stability (|4j), so the flock ultimately follows the leader in equilibrium 
formation. Here we isolate the following problem: What is the transient 
behavior of Zk{t)l The orbit of the agent farthest from the leader, zx{t\ is 
plausibly the one that suffers the worst effects (and this is amply born out by 
numerous numerical experiments). To simplify the discussion we concentrate 
on its orbit. 

This problem constitutes part of a larger research project that aims to 
analytically understand the dynamics of a large number of agents trying 
to move coherently in a changing environment that causes a few elements, 
the leaders, to react to it. An interesting class of systems is given by the 
requirements that most agents "receive information" from the same number 
of nearby agents (homogeneity), and that this information only consists of 
their neighbors' position and velocity. This pattern only changes for agents 
on the boundary of the flock (where it has fewer neighbors). In our case we 
furthermore insist on the interaction being linear and involving only positions 
and velocities of the neighbors. The aim is to give a qualitative analysis 
of the transients of these systems ("Newtonian Networks" is a concept we 
propose hereby) as the number of agents is very large. The analysis below 
reflects some of the complications of this endeavor by examining a (apparently 
simple) paradigm of this idea. 

Notational Conventions: To avoid confusion, we list two important con- 
ventions here. The first is that we assume that both / and g are negative 
reals to insure asymptotic stability (Theorem 12.21) . The second is that we 
define the symbol yfz as the root with angle in the interval [0, it) (branch cut 
along the positive real axis). 

2. Preliminary Results 

In this section we first give the eigenvalues of the matrix M of of Equation 
(CQ). We then give the expression for the frequency response function for 
the trailing agent and discuss its singularities. The following constant will 
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frequently simplify formulae: 

_ 1-p 1 
k = or p = 

P 1 + AC 

In the statement of the next result and that of Proposition 12. 51 we use the 
following equation, where p G (0, 1) and are real variables: 

(2p- l)cot0 = cotiV0 . (6) 

Recall that the matrix P is defined in Equation flS]). 

Proposition 2.1. (f^J) For any p G (0, 1), the matrix P has N distinct 
eigenvalues {Xe}^ 1 : 

i) If p G (0, \\: fori G {0, ...,N- 1}, \ t = 1 - 2y/p{l-p) cos fa, where 

fa e (f > t) so ^ es ®" 

ii) If p G (|, ^^]: Identical to i). 

iii) If p G (^,1): /or^ G {1, . . . , N - 2} ; A, = 1 - 2y/p(l - p) cos fa, 
where fa G (f , S o/, eS ®; A = ^ (i^)""* + (jfrT'') 
and Aat_i = 2 — A . 

One can show (see [H, 0, 0]) that the eigenvalues of M defined in Equation 
([T]) are given by the solutions vi± of 

v 2 - X e gu — Xif — , (7) 

where Xe runs through the spectrum of P. So: 

Theorem 2.2. The eigenvalues of M are 




where Xe runs through the spectrum of P. Because the Xe are contained in 
the interval [0,2] (see Proposition \2. the system is stabilized (or globally 
stable) if and only if both f and g are strictly smaller than zero. 
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Corollary 2.3. The eigenvalues v±£ of M in the complex v plane either lie 
on the circle W+jj | 2 = namely whenever > 1, or else are real numbers 

less than or equal to — j^j. 

Now we turn to the frequency response function of the trailing car when 

Corollary 2.4. (f^j) For p £ (0, 1)\{|} the frequency response function of 
the last agent is given by 

where [i± = /i±(v) = ± (7 ± ^7 2 - 4p(l - pjj and 7 = 7(1/) = f+ f"~f 
As functions of v, the a^iv) in fact are proper rational functions. 



In what follows the location of the roots of the denominator of a at is 
important. Recall that fi + fi_ = n and define the function / : C — > C equal 
to the denominator: 




N-r 

K 



We note that we have two different representations of the set v±i, one in 
terms of eigenvalues and the other in terms of roots \i±£. The locations of 
the zeros fi±£ are described in the following result 

Proposition 2.5. (fjHj) For any positive real number k, the function f has 
2N+2 roots. Two of these are the fixed points of the involution h : fi — >■ \i = - 

and are given by The remaining 2N roots, {pe±} e w S l , have period 2 

under the involution h and are given as follows: 

i) If K > 1: N roots are given by pe+ = v^ e '*i where <pt £ y^, for 
I £ {0, . . . ,N — 1}, solves (T2J); the remaining roots are the images under h 
of these or: y^e"*, respectively. 

ii) If k £ 1): Identical to i). 

iii) If k £ (0,^j): N — 2 roots are given by fie + = ^/He 1 ^ 1 , where (f>£ £ 
f W > ^^n 1 ) , for i E {1, . . . , N - 2}, solves (TJj; iV - 2 are images of these 
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under h; the remaining roots are fiQ + G (a/k, 1) and its images under h and 
multiplication by -1. We have jj 0+ = 1 — ~(1 — k 2 )^" 1 + 0(k 2N ~ 2 ). 



3. Laplace transform and Residues 

In Equation flTJ we set z (t) = 5(t). Set the initial conditions as follows: 
for all k > 1: i%(0_) = Zfc(O) = 0. Then z k (t) equals the Green's function of 
this problem (see [3(): 



The actual impulse response functions z k {t) can be obtained from this by 
twice integrating with the usual initial conditions (ifc(0) = z k (0) = 0). Here 
we calculate the impulse response of the last car (labeled by N). The strategy 
is to perform a residue expansion (or partial fraction expansion) of on{v) 
(given in Corollary 12.41) . 

Considering Corollary 12.41 we write a^iy) = ^£4 as a quotient of poly- 
nomials with degree (p) at least that of degree (q). The zeros of q are the 
eigenvalues of M. Thus according to Theorem 12. 21 the denominator in a^{v) 
has only simple roots located at except when for some i: —4/ = X^g 2 . 
Avoiding that case for simplicity, we have: 



With this proviso, we will calculate all the residues Res(a7v(^), v ±i ), v_i = vi 
and 



The indexing has been chosen so that the pair v±i correspondings to a pair 
\i±i of zeroes of /. 




(8) 





(9) 
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This representation allows us to 'compute' the motion of the N-th agent 
via the inverse Laplace transform 

z N (t) = £ f teMv),vt-) eVl _ t + ^(a N (u),u e+ ) eUl+t \ +Cn + 
i=o \ £- e+ / 

(10) 

The constants of integration Cn and Dm have to guarantee that zn{0) = 0, 
z N {0) = 0. 

Theorem 3.1. If the poles are simple then 

(f + gu i± y ^^{^-Kf 



Res(ajv(V), V£±) 



v l± (2f + gu t± ) 2N^~ 2 ^ 2 - 1) + 2n™ + 2k n ^ 



Proof: If in the above Proposition we replace /i_ by k//j,+, then the expres- 
sion for in Corollary 12.41 is a rational function of p + alone: 



a N 



1 + K K N /J,+ (fl+ — K.) 1 + K Pn(^- 



(The polynomials pn and q7v still have a factor — k) in common, which is 
kept to simplify the calculation.) Recall that p + is a function of 7 by (choose 
the "+"root): 

(1 - p) - 7/i + p/i 2 = , 
with 7 determined by ve± through 

v 2 

7 = 1- 



f + gv 

The pole expansion of is performed as in jjij]. When the poles of a N are 
simple we obtain that 

Res(a N (v),v e ±) = —J — -Res(a N (n), ni±) 

Using the above relations, one obtains: 

/i+Ott) 2 vt±(2f + gut±) 
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Using this and replacing the residue of a^(p + ) by yj^j- 

.^{tili-K} (f + gv e± )' 

Res(a N (is),is e ±) = + 



fJ-Xe J »i± (2/ + gvt±) ' 



2N^~ 1 if^ - 1) + 2^(jm™ + k»-i] 
which after some simplification gives the desired result. 



4. The Dominant Poles 

We show that there are three cases: 

1. When p > 1/2 (or k < 1), the weight is more on the agent following. 
In this case two poles dominate the frequency response a^{v) and we 
can estimate the impulse response by the inverse Laplace transform. 

2. When p < 1/2 (or k > 1), the weight favors the agent in front. In 
this case no poles appear to be negligible, and the inverse transform is 
problematic. 

3. When p = 1/2 (or k — 1), equal weight is on the front and back 
neighbor. In this case on the order \/N poles dominate. The inverse 
transform can be done and this case is described in 

Theorem 12 . 2 1 implies that the sign of 1 + ^ determines whether the eigen- 
values of the system are real or complex. Proposition 12.11 tells us that when 
k < 1 the eigenvalue Ao is exponentially small. Given our assumptions (all 
parameters fixed, except N — > oo) we may thus assume that vq± are complex 
when k < 1. When k > 1 it is certainly possible that all eigenvalues are real. 
To simplify the analysis we will assume from now on that all eigenvalues are 
pairs of complex conjugates (or that 2™ > 1). 

We now present the magnitudes of the residues as well as the relative 
location of the eigenvalues v±e to the real axis in terms of a table. First, 
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define the following factors: 

k n 1 ^j rt ' (a* 2 ^ - *0 2 

2N^f l -\t? +l -l)+2 t i\ N t +2K N -^ 
1 

The factors / and II multiply to produce the residue, while factor III 
describes the relative inverse location of the poles to the imaginary axis: /// 
is large when the poles are near the imaginary axis. The following tables 
hold: 



K < 1 














# 


A 


V 


M 


i 


II 


III 


0+ 






1- l(l_ K 2 )K iV-l 






2 

SoTsl 


0- 




^Aos - «\AoT71 


k(1 + i(l - K 2 )^™- 1 ) 


2%Ao 


A 


2 

A 0l9 


1+ 


1 - 2^/p(l - p) cos 4> t 


O(l) 


0(1) 


O(l) 




O(l) 


I— 


1 - - p) cos <j> t 


O(l) 


O(l) 


O(l) 


0( K N ' 2 ) 


O(l) 


(JV - 1) + 


2 - A 


O(l) 


O(l) 


O(l) 




O(l) 


(JV-1)- 


2 - \o 


O(l) 


O(l) 


0(1) 


0(k jv ) 


O(l) 



ft > 1 














# 


A 


V 




I 


II 


III 




1 - 2\/ p(l — p) cos 


X 1 




X 1 




X 1 



(The symbol / x 1 means that the absolute value of the expression "I" 
is 0(1) but not o(l)). 

The accuracy in the table is as follows. Where the entry equals ax the ac- 
curacy is 0(x 2 ), where the term is ax + bx 2 (such as in both the //-column and 
the I-column) the accuracy is 0(x 3 ). Otherwise exceptions are mentioned. 
For example, one of these, the entry for i>£±, follows because the absolute 
value of the imaginary part of vi± is 0(1). Note that the A's necessarily have 
the same "+" and "-" entries. 

We first discuss the situation when k < 1. For v close enough to the 
pole at is 0+ , a N (u) can be estimated by evaluating the behavior of a N near 
that pole (see Equation 0). Substitute v = i\/Xo\f\ and (see Equation ([9])) 
multiply the factors I, II, and III, of the first line in the table to obtain the 
following result (see also Theorem 4.6 of Q): 

Corollary 4.1. When re < I, a N {i^J\) = -i^- + 0(1). 



I = 

II = 

III = 
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This result should be compared with Theorem 4.6 of [4| (especially its 
proof). 

From Equation fflUj) and the table above, we see that there are only two 
poles that yield an exponential contribution to the impulse response function 
zjv(t). From the table one can calculate this contribution by multiplying I 
and II, and dividing the result by the square of v (see Equation ffHJj) ). 

Corollary 4.2. Using only the two principal poles (when n < I) as an ap- 
proximation, the impulse response for the trailing car is given by: 

z N (t)=t J. e A °^ 2 sin(v^^t) where A = ^(I-k 2 )^ 1 . 

y\j\ vAo 1 



(The two leading poles determine the dynamics of the last agent.) 

The dynamics in this case is virtually entirely determined by the leading 
eigenvalue A of the (reduced) Laplacian. Note that A is exponentially small, 
yet positive, in N. Also notice that the constant g is negative. The motion 
Z]y(t) for a substantial time interval (as long as \/Xo\f \ t is small) is roughly 
equal to ZN(t) = t(l — e Xo9t ^ 2 ) and therefore remains small, that is: for an 
amount of time 0(k~) the agent appears to not move. 

Now we turn to the case where k > 1. From the tables we conclude: 

Corollary 4.3. When n > 1 there are O(N) poles that play a role in the 
dynamics of the orbit of the trailing car. (All have exponentially large residues 
while other factors are 0(1). ) 

The fact that there are no dominant poles in this situation effectively 
prevents us from giving an approximation of the impulse response function. 

5. Impulse Stability 

In [H we suggest the following notion of impulse stability. 

Definition 5.1. Consider Equation with forcing determined by Zo(t) = 
5(t) and subject to the initial conditions Zk(0) = ifc(O) = 0. Let = 
sup t>0 |-C-(z/v(0 — ^o(^)) I ■ The system is called 'impulse stable' if it is asymp- 

d) 1/N 

totically stable and if for i equal to 0, 1 and 2: lim sup^^^ Z N 
Otherwise the system is 'impulse unstable'. 
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Impulse instability in this sense means that if we give the leader a 'unit- 
kick', then that perturbation travels through the flock and causes sup |zjv(0I> 
sup |ijv(t)|, or sup |z/v(t)| to grow exponentially in N, before eventually dying 
out. 

Conjecture 5.2. The system of Equation (TJP is impulse stable if and only 
ifP = \ (or k = 1). 

The fact that for p — | the system is stable in this sense was proved in 

[sj]. When p > |, Corollary 14.21 implies that for an amount of time 0(rz%) 
the trailing agent appears to not move. During all that time however the 
leader has traveled at unit speed. Thus at this point t is in time z^{t) — Zo(t) 
is exponentially large in N. This proves impulse instability when p > ~. 

The problem resides in the case p < |. Due to Corollary 14.31 we cannot 
easily find approximate solutions. Here is a heuristic argument in that case. 
Use the fact that P-P+ = k to rewrite 

M>0 = h±Ji^ ( i - ( a=) ") (id 

K p+-p + \ v+-p+ Kp+J J 

Now one shows (see the appendix of [H) that there is an r e (0, 1) for 
which < r and furthermore that there is an interval (0,a;+) on which 

\p,_{iu)\ > 1. From this it follows that grows exponentially large (in 
N) as N tends to infinity on a fixed interval (0, Thus the norm 

| ajsiiiuS) || 2 of on grows exponentially. Since the Fourier transform (and 
its inverse) preserve the L2 norm (by Plancherel's theorem), we now have 
that zn has an exponentially growing L2 norm. Until this point there is no 
problem. But now we want to prove that the sup- norm of 2jv(t) must also 
grow exponentially. We know that is a linear combination of eigensolu- 
tions each of which decays with e Re ^ Ue± ^. In this case the numbers Re(vi±) 
are uniformly (in N) bounded by a strictly negative number. So it seems 
reasonable that if such a function is to have an exponentially large L 2 norm 
then its sup-norm must grow exponentially as well. However, we have been 
unable to prove this rigorously. 
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